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Abstract 

Total photonuclear absorption cross sections of and ^He are studied using 
realistic NN and NNN forces. Final state interactions are fully included. Two 
NN potential models, the AV14 and the r-space Bonn-A potentials, are con- 
sidered. For the NNN forces the Urbana-VIII and Tucson-Melbourne models 
are employed. We find the cross section to be sensitive to nuclear dynamics. 
Of particular interest in this work is the effect which NNN forces have on the 
cross section. The addition of NNN forces not only lowers the peak height but 
increases the cross section beyond 70 MeV by roughly 15%. Cross sections 
are computed using the Lorentz integral transform method. 
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Typeset using REVT^ 



The process of nuclear photoabsorption is known to be strongly dependent on exchange 
currents and hence on the underlying nuclear dynamics. One therefore expects the photoab- 
sorption cross section to be sensitive to NN and NNN nuclear forces. It is the purpose of 
this paper to explore this sensitivity by computing the total photoabsorption cross sections 
of the trinucleons by using realistic NN and NNN potentials. The only other calculations of 
trinucleon photodisintegraton covering a larger range of energies (0 - 40 MeV) and employing 
realistic NN forces are, to our knowledge, those of Sandhas a/ |T]]. However their calcu- 
lations are for the exclusive process of two-body breakup while we compute the inclusive 
(two-body -|- three-body breakup) cross sections over a larger energy range. It is because we 
are looking for effects of NN and NNN forces that energies in the tail region, >70 MeV, 
are included here. 

We use the Argonne AV14 potential and the Bonn-A r-space potential ||^. As NNN 
potentials we include the Urbana-VIII (UrbVIII) 0] and Tucson-Melbourne (TM) [g NNN 
models. The total photoabsorption cross section for energies below pion threshold is com- 
puted by using the unretarded dipole operator 

z 

D = Y^{f,- R,m) . (1) 
i=l 

This form of the operator already includes the effects of exchange currents due to Siegert's 
theorem. It is known that as far as the total cross section is concerned corrections to the 
unretarded dipole operator are small in the energy range studied here. We consider this 
approximation adequate for investigating the effects of NN and NNN forces to the process. 
The total photoabsorption cross section is then given in terms of the dipole response function 
Rhj 

atot = ATi\e^/nc)E^R{E^) (2) 

where 

R{uj) = < ^^\DJ{H -Eo- uj)D,\^o > ■ (3) 

To calculate R{uj) we use the Lorentz integral transform method introduced in This tech- 
nique eliminates the need to compute final-state continuum wave functions. The transform, 
$((T/j,cr/), of R{uj) is defined by 

H^R,^i) = r du- ^^"^ , . (4) 

It can be expressed as the norm of the square integrable function \1/ 

= < > (5) 

where 

{H - E^^aR + iai)\^ > = D,\^Q> . (6) 

Finally R{ijj) is obtained by inverting the transform (4). The inversion procedure is elab- 
orated upon in This method was applied to the study of few-nucleon responses in 
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||8HT2|. In Faddeev techniques were applied to solve equations of the form of Eq.(6) 



whereas H T3"12| used expansions in a set of correlated hyperspherical harmonics (CHH) 



for this purpose. The latter technique is also adopted in the present work. We expand the 
functions, \l/o and \E', in a set of correlated hyperspherical harmonics (CHH) according to 

= ^ Q 0, (7) 

where a) is a correlation operator and (pi are a totally antisymmetric basis set constructed 
from a spatial part Xi,n ^^id a spin-isospin part O^j^ 

Ct>^ = T. XiJ, ■ (8) 

Whereas in [^,[TD|,|TT| , where simpler potential models were employed, it was sufficient to 
take uj in the state independent form 

= n /(^^.) • (9) 

we now use a state dependent correlation operator of the form 

^ = SX{Y.fsT{n,)PsT{ij) (10) 

i<i S,T 

where Psxiij) are projection operators onto nucleon pairs (ij) with spin S and isospin T and 
where iS is a particle symmetrization operator. This form of correlation operator is easily 
incorporated into the calculations by first constructing the set of coefficients < 9^j_i\uj\9^ > 
and then writing 

= E X^,,' 0,' (11) 



where 



Xi,ti' = E < ^M'l^l^M > Xi,ii ■ (12) 



The correlation functions fsT{f) are chosen as follows. For r < tq , the healing distance, 
fsri^) is chosen to be the zero energy pair wave function in the corresponding ST state. 
Healing is insured by imposing the conditions fsT{f)=^ for r > tq and f'sTiTo)—^- The 
ST= 13 and 31 cases are determined from the ^Sq and ^Si partial waves of the NN potential. 
However for the S'T=11 and 33 cases the ^Pi and '^Pi potentials are not sufficiently attractive 
to obtain a healing distance. Therefore we introduce an additional intermediate range central 
interaction such that a healing distance of 10 fm is obtained. Further details can be found 
in ig. 

Table 1 gives our results for the ground state properties of '^H using the various potential 
models described above. We point out that our results for both the ground state and the 
Lorentz integral transform arise from fully converged expansions in the CHH basis. The 
inversion leads to stable results for R{uj). In order to obtain the correct binding energy 
with the TM NNN potential we have adjusted the cut-off mass A in the monopole form 
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factor. This requires A= 4.67/i and A=4.07/i for use with the AV14 and Bonn-A (r-space) 
cases respectively. Previously published ground state properties are available [Il3i|l4| for the 
AV14 and AV14+UrbVIII potentials. The corresponding results in Table 1 are in agreement 
with these. Since the method of the Lorentz integral transform allows one to calculate the 
cross section over the entire energy range, sum rules can be used to give additional checks. 
An example is provided by explicit evaluation of the inverse energy weighted sum rule 
(T_i = E~^(Ttot{E-y)dE.y which is related to the triton point proton radius through 

This sum rule is calculated in two ways: (i) by direct integration of the response giving 
cr!!\*, and (ii) by use of the ground state wave function giving a^Ji. As seen in Table 1 
agreement is obtained to better than 0.5%. This level of agreement helps confirm that the 
lower energy regions, and especially the peak heights, are correct. A check on the higher 
energy regions is provided by evaluating the TRK sum rule, again by direct integration and 
by the ground state expectation value of the appropriate double commutator. This was 
done for the AV14 model, yielding 70.7 MeV-mb by direct integration and 71.6 MeV-mb by 
evaluating the double commutator. The direct integration was carried out to 2000 MeV. 
It should be noted that integrating up to 300 MeV excitation only exhausts about 90% of 
the sum rule. This is in contrast to the integration of the TRK sum rule produced by a 



soft-core potential such as TRSB potential [|T5| where almost all the strength is exhausted 
at 300 MeV excitation §. 

The value of the cross section at the peak depends on ground state details, especially 
the radius. This is evident from the remarks above on the inverse energy weighted sum 



rule. Fig.l, for example, shows results obtained previously ||T0| using the Malfliet-Tjon 



(MT) |jT6| and Trento (TN) [§] potentials together with new results obtained here for the 
AV14 and AV14+UrbVIII potentials . Ground state point proton radii computed from these 
potentials are 1.61, 1.59, 1.66, and 1.60 fm respectively. The similar peak heights for these 
rather different potentials is due largely to the similarity of these radii. However, we have 
observed that in all cases considered here the combination of realistic NN forces with NNN 
forces leads to roughly a 10 % decrease in the peak height. Moreover, even though the four 
potentials give similar radii and binding energies the AV14 and AV14+UrbVIII potentials 
produce cross sections with significantly larger high energy tails than the semi-realistic MT 
and TN models. Part of the difference must be attributable to the tensor force contained in 
the realistic potential models. 

Figs. 2(a) and 2(b) show the peak and tail regions of the T=l/2 contribution to the 
absorption cross section as computed from the five potential models of Table 1. A study of 
these figures reveals that the addition of three-body forces (3BF) lowers the peak height and 
increases the cross section in the tail region by 10 - 20%. The lowering of the peak can be 
understood since the addition of 3BF increases the binding energy and decreases the radius 
in all cases. Despite the increase of the cross section in the tail region due to 3BF, no clear 
separation between the cross section with and without 3BF is seen in this channel. We note 
the importance of binding in this T=l/2 channel by observing in the tail region the rather 
disparate curves for the pure NN cases of AV14 and Bonn-A. Recall from Table 1 that these 
NN potentials correspond to binding energies of 7.69 MeV and 8.15 MeV respectively. In 
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addition the three curves corresponding to the addition of NNN forces are rather close and 
difficult to distinguish, reinforcing the argument that effects here are mainly due to binding. 

Curves corresponding to Fig. 2 but for the T=3/2 channel are shown in Figs. 3(a) and 
3(b). The behaviour of the cross section in the peak region is similar to the T=l/2 case 
although here the peak height is less sensitive to the interaction model. Fig. 3(b) shows that 
the T=3/2 channel clearly separates all models containing a 3BF from those with purely NN 
forces. This is a 5 - 10 % effect in the energy range 70-110 MeV. However contrary to the 
T=l/2 case, binding effects here do not play an important role as evidenced by the nearly 
overlapping NN curves in the tail region. As a result the separation between the curves in 
Fig. 3(b) is dominantly an effect due to NNN forces i.e. the T=3/2 channel is more sensitive 
to NNN forces in the tail region. The import of this observation is that it is the three-body 
breakup channels which might show more evidence of 3BF since they carry all of the T=3/2 
strength. In fact analyzing our theoretical results and the experimental nd and ppn data 
from Ref. |TS| leads us to conclude that in this energy region the T=l/2 content of the ppn 
channel is about 1/3. Our findings are in line with proposals Jl^ which suggest particular 



kinematical set-ups in the nnp or npp channels to be sensitive to the presence of 3BF. 

Figs. 4(a) and 4(b) compare to available data p^p!9| the total cross sections for and 
■^He as computed with two of our more complete potential models, AV14+UrbVIII and Bonn- 
A + TM. The '^He calculations include the full Coulomb interaction between the protons. 
Our calculations agree rather closely with the available data. We note by reference to Fig. 
1 that the semi-realistic models do not yield the same level of agreement with the data. 

We have computed the total photo-absorption cross sections of the trinucleons by using 
realistic NN and NNN forces. Final state interactions have been fully included through 
the use of the Lorentz integral transform method. The results show clear sensitivity to the 
underlying nuclear dynamics. In particular the tail region of the T=3/2 channel appears 
to be more sensitive to details of the NNN forces than to NN forces. Therefore further 
theoretical and experimental effort should be devoted to this promising kinematical region. 
On the theoretical side several corrections, expected to be minor, will be investigated in the 
future. These are the effects of higher multipoles, retardation in the one-body operators 
and inclusion of explicit two-body currents beyond the Siegert theorem. In addition the 



calculation will be extended to the new charge dependent potentials such as the AV18 ||20 
and Nijmegen |21[] potentials together with corresponding NNN models. The rather large 
effect of the NNN force in the tail region could imply a much larger effect in some selected 
kinematics of the exclusive three-body breakup reaction. In this regard it should be pointed 
out that a modification of the Lorentz integral transform method used here p^ , |23[| will allow 
the calculation of exclusive reaction cross sections as well. 
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TABLES 

TABLE I. Ground state properties for AV14, Bonn-A, AV14+UrbVIII, AV14+TM, and 
Bonn-A+TM potential models. 



Potential 


Eb (MeV) 


yj< rl > (fm) 


Pd 


ai"/ (mb) 


a^i (mb) 


AV14 


7.69 


1.66 


8.94 


2.636 


2.645 


Bonn-A 


8.15 


1.61 


6.88 


2.481 


2.486 


AV14+UrbVIII 


8.49 


1.60 


9.67 


2.442 


2.453 


AV14+TM 


8.48 


1.60 


9.67 


2.448 


2.459 


Bonn-A +TM 


8.47 


1.59 


7.25 


2.411 


2.415 
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FIGURES 



FIG. 1. Total photoabsorption cross sections predicted by the TN (dotted curve), MT 
(dashed curve), AV14 (dash-dot curve) and the AV14 + UrbVIII (solid curve) potential models. 

FIG. 2. T=l/2 part of the total photoabsorption cross sections in the peak region (a) and 
tail region (b). The corresponding potentials are AV14 (dotted), AV14 + UrbVIII (short dash), 
AV14 + TM (long dash), r-space Bonn- A (dash-dot), r-space Bonn-A -|- TM (solid). 

FIG. 3. As in Fig. 2 but for T=3/2. 

FIG. 4. Total photoabsorption cross sections for "^H (a) and '^Hc (b) computed from the AV14 
+ UrbVIII (solid) and r-space Bonn-A + TM (dashed) models. The shaded area represents the 
data from [19] while the black dots are the data of [18]. 
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Fig. 4 



